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ABSTRACT 

The 4050A band of C3 was observed with Keck/HIRES echelle spectrometer during 
the Deep Impact encounter. We perform a 2-dimensional analysis of the exposures in 
order to study the spatial, spectral, and temporal changes in the emission spectrum 
of C3. The rotational population distribution changes after impact, beginning with an 
excitation temperature of ~45 K at impact and increasing for 2 hr up to a maximum of 
61±5 K. From 2 to 4 hours after impact, the excitation temperature decreases to the pre- 
impact value. We measured the quiescent production rate of C3 before the encounter 
to be 1.0 x 10 23 s _1 , while 2 hours after impact we recorded a peak production rate of 
1.7 x 10 23 s _1 . Whereas the excitation temperature returned to the pre-impact value 
during the observations, the production rate remained elevated, decreasing slowly, until 
the end of the 4hr observations. These results are interpreted in terms of changing 
gas densities in the coma and short-term changes in the primary chemical production 
mechanism for C3. 

Subject headings: comets: general — comets: individual (9P/Tempel 1) — molecular 
processes 



1. Introduction 



The Deep Impact (DI) mission created an artificial outburst on comet 9P/Tempel 1 by exca- 



vating material from the subsurface (A'Hearn et al. 2005). The event was observed by a network 



of orbiting (e.g., Feldman et al. (2006); Lisse et al. (2006)) and ground-based facilities (Meech 



et al. 2005). iFeldman et al. (2006) use ultraviolet observations to show that the amount of CO 



relative to H2O is similar in both the impact ejecta and in the outgassing of the quiescent comet. 



Jehin et al. (2006) measure 12 C/ 13 C and 14 N/ 15 N in the material released by impact and find it to 



be the same as the surface material. Mumma et al. (2005) show that the relative abundances of 



CH3OH and HCN are the same before and after impact, although the amount of C2H6 relative to 

a factor of two shortly after impact 
including the production rate of C3 



H2O increases by almost a factor of two shortly after impact. Schleicher et al. (]2006|) report that 
in nearly all respects 



comet 9P/Tempel 1 returned to 
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pre-impact conditions within 6 days. With the precisely timed outburst caused by DI, the short 
(4 hr) time period directly after impact offers a unique opportunity to observe molecular formation 
mechanisms that occur on short timescales. Indeed, Jackson et al. (2009) use two- or three-step 
Haser models to study temporal changes in the total emission from the following species: O, OH, 
CN, C2, C3, NH, and NH2. Here we describe the observed variations in the excitation temperature 
and production rate of C3 in the first four hours after DI. 



Observations and Reduction 



The observations were conducted using the HIRES cross-dispersed echelle spectrometer (Vogt 



et al. 1994) at the W.M. Keck observatory. HIRES offers nearly complete spectral coverage from 
3000 - 5880 A at a spectral resolution, A/AA ~ 48,000 with a 7'.'0 x / . / 86 slit. The plate scale is 
/ /239. Three 20min exposures were taken on 30 May 2005 UT and combined to characterize the 
comet before impact. On 04 Jul 2005 UT, a series of shorter exposures were taken to capture short 
timescale changes in the spectrum after impact. Exposure times increase with airmass and range 
from 10 - 30 min during the first 3 hrs after impact. An analysis and description of the observations 
is given in Jehin et al. (2006), with additional details in Cochran et al.l (2007) and Jackson et al. 



(2009). Here we focus on the 4050 A band of C3 (from 4049 - 4057 A), which falls onto the central 
1/8 of an echelle order, near the peak of aperture blaze. 

Standard procedures are used to bias correct and flat field each target and calibration exposure. 
Cosmic rays are removed using a two-step process. First, large cosmic ray hits are removed by 
considering the 13 target exposures together at each pixel and replacing large deviations from the 
mean with the median value. Then, to identify smaller cosmic ray hits, a 6x6 pixel spatial filter 
is used on each image, where pixels exceeding the mean value within the box by 2.5cr are replaced 
by the mean value. The echelle order is linear along the CCD over this wavelength range and it is 
sufficient to rectify the spectra by rotating the image 1.15°. Rectification is verified to be better 
than a pixel by comparing the dust continuum spatial profiles along the order. Flux calibration is 
performed for both the pre-impact and post-impact exposures using exposures of the flux standards, 



Feige 67 and BD+284211, respectively (Oke 1990). The total counts are integrated along the slit 



for the calibration standard and averaged over the bandpass of the Oke ( 1990 ) observations (4020- 
4060 A) . The HIRES observations confirm that there are no significant spectral features in either of 
the calibration stars for this bandpass, so the mean observed count rates and the literature values 
for the total flux at 4050 A are used to convert from observed count rate (DN /s) to flux density (ergs 
s cm~ 2 Hz -1 ). We correct for the different airmass during the target and calibration exposures 
using a characteristic atmospheric extinction, k = 0.56, at 4000 A. 

The solar continuum dominates the cometary spectrum at most of the sampled wavelengths, 
with temporal variations in continuum reflectivity due to the changing particle size and composi- 
tion of the ejected material. We use the median value along the dispersion axis to calculate the 
spatial profile of reflected sunlight at each time step to empirically remove the solar contribution. 
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The median profile is convolved with a high-resolution, R ~ 85, 000, UVES solar spectrum^] and 
interpolated onto the HIRES plate scale to model the reflected sunlight. In the narrow wavelength 
range being considered, the changes in the reflectivity spectrum of the ejecta are small and all 
observations are successfully fit with one spectral template. This model is subtracted from the 
observations, leaving only the gas emission. The higher signal-to-noise spectra taken on 30 May 
2005 UT are used to demonstrate the dust continuum removal, Figure 1. Before removing the solar 
spectrum, some transitions in the emission spectrum of C3 (e.g., near 4052.5 and 4053.2 A) are not 
apparent in the observations because they are collocated with solar features. After removing the 
solar spectrum these lines stand out against the background. 



3. Temperatures and Production Rates 



The C3 excitation profile — that is, the distribution of population among rotational levels - 
can depend on the temperature and density in the coma, the heliocentric velocity and distance, 
and the formation mechanism of the molecule (e.g., the parent and grandparent molecules of C3). 



Rousselot et al. (2001) describe a statistical equilibrium model that can be used to interpret the 



rotationally-resolved emission spectrum of cometary C3 by taking into account changes in helio- 
centric velocities and distances. However, changes in those values are small on the timescale of 
these observations and we assume a simple thermal distribution of states to model the emission 
spectrum. This approach is sufficient for the purposes of measuring the relative changes in the 
rotational excitation profile among exposures. 



Tanabashi et al. 



(2005) for the A x Il u 



Spectra of C3 are calculated using the line list of 
(000-000) transitions. The fraction of the total molecules that are in a particular J level is given 
by 

2J + 1 f-hcB 



Fj(T r ) 



-exp 



where q r is the rotational partition function given by 
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For C3 the ground state rotational constant Bq=0. 43057 cm 1 ( Schmuttenmaer et al. 1990) 



Each modeled spectrum (e.g., Figure 2) has three free parameters; the excitation temperature, 
T r , line width, w, and an intensity scaling factor, a, which are used to reproduce the observed 
spectra (Table 1). Spectra are calculated for a range of T r , w, and a and compared to obser- 
vations. The quality of fit is determined by \ 2 = Ylii^Vi l a2 \ where dyi is the residual at each 
sampled wavelength, i, a is approximated by the noise in the observed spectrum, and the best fit 



2 http:/ /www. eso.org/observing/dfo/quality/UVES/pipeline/solar _spectrum.html 
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corresponds to Xmin- Locations in parameter space where x 2 ~ X m in = 1 t ne uncertainties in 
the input parameters. The residuals from the fit indicate that the thermal population distribution 
is moderately successful at reproducing the observations. A successful fit would be indicated by 
residuals that are random noise, however there is structure in the residuals due to the non-thermal 
population distribution. 

We integrate the flux from C3 over 4049 - 4058 A and use the standard fluorescence efficiency 
factor of g (r#) = 1.0 x 10~ 12 x rjj 2 , at a heliocentric distance, r# = 1.55 AU ( A'Hearn et al.|[l995 ), 
to find the total number of C3 molecules in the sampled area. A Haser model ( Haser|[l957 ) assumes 
isotropic outgassing from the nucleus at a constant velocity, v. Parent species are destroyed to form 
radicals with lifetimes of t p and r r , which define characteristic length scales l p = t p v and l r = r r v 
for the parent and radical, respectively. The parent species and radical are assumed to travel in 
the same direction and at the same speed. The number density of species at a distance, x, from 
the nucleus is given by 



nixj 



Q 



Attvx 2 \ L — I 



(3) 



As detailed in Newburn Sz Spinrad (1984) the relationship between the production rate, Q (s ) 
and the column density profile, N(x) (cm -2 ), is given by 



N(x) 



Q 2 / l r 



Attv x V L — I 



x/l r 



ko(y)dy 



x/l p 



ko(y)dy 



(4) 



where x is the distance from nucleus (cm), and ko is the zero-order modified Bessel function of the 
second kind. 

The projected area that we observe is small, so these observations are not useful for constraining 
scale lengths and we use the values of l p = 6.7 x 10 3 km and l r = 6.5 x 10 4 km at a rfl-=1.55AU 



(A'Hearn et al. 1995). 



The spectra taken were recorded in 13 exposures over the first 4hrs after impact (until the 
comet set), and when extracted along the slit show temporal changes in both excitation profile 
(Figure 3A), and the total flux ( [Cochran et al. 2007; Jackson et al. 2009). The high- J lines in 
the i?-branch band-head (near 4050 A) are more prominent in the exposures roughly 100 min after 
impact (Figure 3A). The spatial profiles along the slit also change shortly after impact (Figure 3B). 
Three examples of the radial dependence of the C3 column density are shown in Figures 3C - 3E, 
along with the best fit production rates, Q, and corresponding model profiles. The C3 production 
rates determined using this technique fall in the range of 1.1 - 1.7xl0 23 s _1 , which are consistent 



with the pre-impact measurements of Lara et al. (2006), however they are an order of magnitude 



smaller than the imaging measurements of Schleicher et al. (2006) and the analysis of Cochran 



et al. (2009). Discrepancies are primarily due to the choice of scale length, l r , which is a factor of 



5.5 smaller than used by Cochran et al. ( |2009 ). The temporal changes in the C3 production rates 
directly after impact are shown in Figure 4 (bottom). 
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The trend in production rates is generally consistent with integrated flux of C3 emission pre- 

In both datasets there are no changes over 
pre-impact conditions for the first 50min (3000s) of observations before a monotonic increase in 
production rates (or integrated flux) over ~130min (7800s). The peak in production rate and in- 
tegrated flux are a factor of ~2.5 larger than the pre-impact values. After the peak there is a very 
gradual decrease until the end of the observations. 



sented in Figure 5 of |Jackson fe Cochran| (2009 



Discussion &; Conclusions 



Since the heliocentric velocity of Tempel 1 does not change significantly over the 4 hrs of the 
observations, the excitation spectrum does not change due to the Swings effect (e.g., changes in 
the radiative excitation of C3 caused by a Doppler shift of absorption lines in the solar spectrum) 
and some other process must account for the observed change in the excitation temperature, T r , 
after impact (Figure 4). Variations in the gas density around the ejecta are one way to explain 
changes in the excitation profile. C3 lacks a permanent dipole moment so radiative relaxation is 
not an efficient method for removing population from high-J states. At high densities however, 
collisional de-excitation can thermalize the population distribution. 



Adamkovics et al. 



(2003) have 



shown that in the interstellar medium the C3 excitation profile depends on density, such that T r 
exceeds the kinetic temperature at densities below 500 cm -3 . 

If the gas density in the coma directly after impact is large, then T r is essentially indicative of 
the thermal temperature. We measured T r ~45K right before impact, which then increased up to 
60 K at 100 min after impact — the same temperature it was in May. This increase may be due to 
decreasing gas density and the lack of collisional de-excitation. However, there is then a puzzling 
decrease in T r from 100 - 200 min after impact. If the gas density were monotonically decreasing, 
then T r should increase and then plateau. However, there is the possibility that the 15 K change 
in T r occurs independently of the changes caused by DI. The quiescent T r measured on 30 May 
2005 is the same as the peak post-impact value of 60 K, so that perhaps T r varies with time for C3 
. The decrease in T r at times greater than 100 min after impact supports the possibility that T r is 
gradually fluctuating on a timescale of hours. 

The production rate for C3 reaches a plateau 130 min after DI and then appears to decline 
after ~ 183 min. The relative flux of CN follows a similar progression with time (Jehin et al 



2006), however CN reaches maximum production at a mid-exposure time of 96 min, significantly 
earlier than C3. One simplistic interpretation is that there are additional intermediate reactions 
between the photodissociation of the parent molecule and the formation of C3. The primary 
production pathway of cometary C3 is the photodissociation of either propyne (H3CCCH) or allene 
(H2CCCH2) — in either case, one of the isomers of C3H4 leads to the formation of C3 via the C3H2 



radical intermediate (Helbert et al. 20051. This common radical means that multiple pathways 



can produce the same rotational excitation spectrum of C3 (Song et al. 1994) and so the parent 



molecule of C3 cannot be distinguished. Jackson et al. (2009) use a 3-step chemical model with 
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the photo-dissociation of the C3H2 radical as the production pathway for C3, and consider either 
allene or propyne as the precursor to the parent radical. However, other mechanisms have been 
hypothesized for the formation of C3, which could produce C3 with a different excitation spectrum. 



Helbert et al. (2005) mention the electron impact dissociation of C3H4 as a source of C3 but 



note that rates for individual reactions have not been determined and hence the relevance of this 
mechanism is uncertain. The propynal radical (C3H2O) has also been proposed as a source of C3 



( Krasnopolsky 1991), and proceeds via an excited state intermediate, C3H2. This radical could be 
the parent molecule of C3 with a different T r than when produced from C3H2. In general, the total 
yield of C3 from C3H2O is only ~1% of the production from C3H2, yet this mechanism may have 
an increased relevance in the 1-2 hrs after DI. Similarly, the dissociative recombination of C3H;j~ 
may yield C3 with a T r that is larger than when produced by the photodissociation of C3H2. It is 
unclear if there is link between the fact that Q{C%) stops increasing 2 hrs after impact, roughly the 
same time that T r starts to decrease, however these two profiles together provide a constraint on 
the chemistry of carbon-bearing molecules during DI. The excitation profile of C3 thus serves as a 
unique diagnostic of the chemistry and physical processes on short length-scales. 

Future studies could test the effect of variations in the gas density on the excitation temperature 
by using an analogous molecule, such as C2, to make an independent measurement of the gas 
density. The 5100 A Swan bands of C2 are also recorded in the publicly available Keck spectra, 



which together with models such as those of Gredel et al. (1989) and Rousselot et al. (2001) could 



be used to quantitatively compare the excitation profiles of both C2 and C3. Such an analysis would 
also provide a constraint on chemical models of carbon-bearing molecules in comets, since variations 
in the formation mechanism of C3 can change the excitation temperature. Detailed studies of the 
time-dependent kinetics of both C2 and C3 will shed light on mechanisms responsible for variations 
in the excitation profile. 



The data presented in this paper were obtained at the W. M. Keck Observatory, which was 
made possible by the financial support of the W. M. Keck Foundation. The authors wish to 
recognize and acknowledge the very significant cultural role and reverence that the summit of 
Mauna Kea has always had with the indigenous Hawaiian community. 
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Time 




w 


a 


8.3 


48.2 ± 7.4 


0.072 ± 0.015 


1.02 


19.2 


44.8 ± 8.8 


0.066 ± 0.021 


0.72 


32.6 


46.2 ± 6.9 


0.059 ± 0.012 


0.76 


48.4 


47.0 ± 4.9 


0.051 ± 0.006 


0.88 


64.3 


48.1 ± 4.2 


0.054 ± 0.005 


0.87 


80.2 


53.7 ± 4.1 


0.062 ± 0.006 


0.83 


96.2 


60.4 ± 4.5 


0.067 ± 0.006 


0.84 


112.1 


55.3 ± 3.6 


0.065 ± 0.005 


0.92 


128.0 


54.0 ± 3.9 


0.062 ± 0.006 


0.89 


143.9 


50.8 ± 3.8 


0.064 ± 0.006 


0.81 


159.8 


46.5 ± 3.8 


0.063 ± 0.006 


0.87 


183.2 


46.9 ± 4.2 


0.064 ± 0.007 


0.85 


214.1 


45.1 ± 5.1 


0.065 ± 0.009 


0.85 



Table 1 : Fit parameters for all spectra, which are identified by time from impact to mid-exposure in 
minutes. Excitation temperatures T r are presented in K, line widths w in A, as well as the unitless 
scaling parameter, a. 



This preprint was prepared with the AAS IfTjrjX macros v5.2. 
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Fig. 1. — An example of the 2-dimensional solar continuum subtraction using the 30 May 2005 
UT spectrum of comet 9P/Tempel 1 (obs.). The y-axis is the spatial dimension along the 7" spec- 
trometer slit. The modeled solar spectrum (see text) is in the middle (model) and the continuum- 
subtracted 2-d emission spectrum of C3 is shown at the bottom (C3). 
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Fig. 2. — The spectrum of C3 before DI (on 30 May 2005 UT) integrated along the slit, with a best 
fit thermal excitation (T r = 60±3K) model (blue) for the rotational population distribution (top). 
The observed spectrum is the average of three 20 min exposures. Residuals between the observations 
and the model, along with laboratory assignments of the transitions from two long-lived perturbing 
states (bottom). 
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Fig. 3. — (A) Time series of spectra extracted along the slit. Spectra are labeled by the time 
at mid-exposure from DI in minutes. (B) The 2-dimensional spectra of C3 after modeling and 
subtracting the contribution of reflected sunlight off dust. (C)-(E) Examples of the column density 
of C3, N(x), as a function of the distance from the nucleus (black lines). For reference, a column 
density profile calculated using the pre-impact production rate of C3 (Q=10 23 mol/s) is plotted 
with a dashed red curve in each panel. The best fit profile calculated using a Haser model with a 
variable Q(Cs) is plotted with a solid blue curve. 
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Fig. 4. — (Upper) Comparison of the changes in C3 excitation temperature and production rate 
after impact. T r is measured by fitting a thermal population distribution to the observed spectra 
and is plotted here as a function of the mid-exposure time after impact. See text for calculation 
of uncertainties. Error bars along the x-axis indicate the length of the exposure. (Lower) The 
production rate, Q(Cs), after DI determined by fitting the C3 column density profile, N(x), with a 
Haser model (see Figure 3 for profiles). Left arrows indicate pre-impact values from 30 May 2005 
UT. 



